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We study the competition between different possible ground states of the double-exchange model 
with strong ferromagnetic exchange interaction between itinerant electrons and local spins. Both 
for classical and quantum treatment of the local spins the homogeneous canted state is shown to be 
unstable against a phase separation. The conditions for the phase separation into the mixture of 
the antiferromagnetic and ferromagnetic/canted states are given. We also discuss another possible 
realization of the phase-separated state: ferromagnetic polarons embedded into an antiferromagnetic 
surrounding. The general picture of a percolated state, which emerges from these considerations, is 
discussed and compared with results of recent experiments on doped manganaties. 



PACS numbers: 75.10.-b, 75.30.Mb, 7530.Kz 



I. INTRODUCTION 

The double-exchange model describing itinerant electrons interacting with local spins was first used to explain 
ferromagnetism in metals, such as Ni and Fe pL The revival of interest to this model was prompted by the recent 
observation of the colossal magnetoresistance (CMR) in manganites The limit of strong ferromagnetic (FM) 
exchange between the conducting electrons and spins, relevant for manganites, was first discussed by de Gennes ||. 
He suggested that the competition between the antiferromagnetic (AFM) superexchange and the double exchange 
results in the canting of the AFM state (the angle 8 between the spins from different sublattices becomes smaller than 
7r). The canting angle grows with the concentration of charge carriers, which explains the increase of magnetization 
upon doping observed in Lai-^Ca^MnOs. 

However, already rather long ago arguments against the canted ground state were put forward, most notably by 
Nagaev [|| . In the de Gennes approach the local spins were treated classically. Quantum corrections stabilize the AF 
state and the canting appears only above certain concentration of charge carriers [0. 

A more fundamental problem is that at partial filling of the conduction band any homogeneous ground state may 
be unstable against a phase separation. Such an instability takes place in the Hubbard model away from half-filling 
Pj and in the t — J model JH|. In a wide range of parameters the ground state of these models is phase-separated, 
stable homogeneous solutions being rather an exception. In this respect the double-exchange model is rather similar, 
which can be seen already from the tendency to a formation of magnetic polarons in this model . These polarons 
are FM droplets surrounding charge carriers in the AF background. If for some reason (large effective mass, disorder) 
the polarons become immobile, they can be viewed as a form of the FM-AFM phase separation (see also [|]j|]). 

The phase separation in the double-exchange model would imply that many experimental data on doped mangan- 
ites should be reinterpreted taking into acount the magnetic and electronic inhomogcneity of the ground state. In 
particular, the charge transport and the metal-insulator transition should be described in terms of percolation rather 
than by the properties of some pure states. The percolation in managanites was emphasized by Gorkov and Kresin 
in Ref. fic(| . Some recent experimental results also strongly point in this direction [pd|-^3[. 

In this paper we discuss the instability against the phase separation in the double-exchange model. Although we 
mostly keep in mind applications to CMR manganites, our treatment has a general character. We go beyond the 
classical treatment of local spins used in most numerical studies of the double-exchange model. Our treatment has 
some points in common with the previous investigations. Nevertheless, we think it is worthwhile to include (with 
necessary modifications) some old results, so as to have a full coherent picture of the behavior of this system. 

Our paper is organized as follows. After formulating the model in Sec. 0, we discuss the homogeneous canted 
state, taking into account quantum effects. Several critical concentrations separating different regimes are identified. 



In Sec. Ill we demonstrate that for small concentrations of charge carriers the homogeneous canted state is unstable 



against a phase separation. In Sec. IV the general phase diagram of the double-exchange model is discussed. The 
Maxwell construction is used to investigate the characteristics of the coexisting phases. It is shown, in particular, 
that depending on the parameters, there may be phase separation not only on AFM and FM phases, but also on 
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AFM and canted states. Then in Sec. [v|, we consider the polaronic state and explicitely show that the energy of this 
inhomogeneous state is lower than the energy of the homogeneous canted state, even when quantum effects are taken 
into account. Finally, in the concluding Sec. VI we summarize the main results and discuss the resulting picture in 
relation to some recent experimental observations. 



II. ENERGY OF HOMOGENEOUS CANTED STATE 

We consider the double-exchange Hamiltonian describing the interaction of itinerant electrons with localized spins, 
to which we add the Heisenberg interaction between the spins (FM Kondo-lattice model): 
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Here, the first term describes the hopping of conduction electrons, the second term is the on-site FM exchange between 
the spin of a conduction electron s and the localized spin S, and the last term is the AFM exchange between the 
nearest localized spins (Jh, J > 0). We assume the value of localized spins to be large: S > 1. Futhermore, we 
consider the case of strong Hund coupling Jh and weak Heisenberg coupling J: 



J H S > t > JS 2 . 



(2) 



Finally, in this paper we consider the simple cubic lattice. 

The model Eq.(Q) describes the competition between the direct antiferromagnetic exchange of the localized spins 
and the double exchange via the conduction electrons, which tends to order the spins ferromagnetically. The simplest 
state that could result from this competition, is a canted state, which has the features of both the FM and AFM 
states: as in the AFM state, the localized spins in the canted state form two sublattices, however an angle 8 between 
the magnetizations of the sublattices is, in general, arbitrary, so that the total magnetic momentum of the system is 
nonzero. The canted state interpolates between the AFM state (0 = tt) and the FM state (9 = 0). 

In the classical treatment of localized spins, used in Ref. ra, the orientation of the spins is fixed: 
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where A and B denote the two sublattices and the unit vectors Ha and ub describe the direction of the magnetization 
of the corresponding sublattice. Such a state we shall call the classical canted (CC) state. 

In another approach Q , one assumes that Eq. (^) only holds on sites that are not occupied by conduction electrons. 
In the quantum language, the localized spins on empty sites have the maximal projection, +S, on the magnetization 
vector of the corresponding sublattice. On occupied sites, however, the localized spin S and the spin of a conduction 
electron s form a state with the maximal total spin Stot = S + ^, because it has the lowest energy in the large- Jh limit. 
The projection of the total spin on the axis of the corresponding sublattice can have only two values: Sf ot = S ± | 
(otherwise, the localized spin cannot remain in the state with the maximal projection after the conduction electron 
leaves the site) . The hopping of electrons between the two sublattices is then described by a 2 x 2 matrix. Diagonalizing 
this matrix one obtains two electron bands with the corresponding hopping amplitudes B, 
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where 7 is defined by sinh7 = S cos %j\j2S + 1. The two-band canted state is refered below as the quantum canted 
(QC) state (though the quantum fluctuations of local spins are not taken into account in this approach). 

The origin of the two hopping amplitudes t± can be easily understood on the simple examples. For the AFM state 
(6 = tt) the two have the same width: t + = t_ = / 2 s+i • ^ n electron, propagating coherently in the AFM spin 
background, creates a state with 5 t z ot = S + h on one sublattice and a state with Sf ot = S 



i on the other: 
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which closely resembles a type of motion considered by Zaanen et al. in Ref.( Similarity, for the FM ordering 

(8 = 0) we obtain from (Q) i + = t and t- — 2 s+\ ■ Here, t + corresponds to the motion of a spin- up electron in the 
spin-up FM background, whereas t_ describes the motion of the spin-down electron in the same background. In the 
latter case, on each site the electron spin and the local spin form the \S^ ot — S — i) state. 



We denote by n± the densities of electrons occupying the two bands, so that 

n+ + n_ — n, 



(5) 



where n is the total electron density. As we shall shortly see, under condition (Q) the transitions between the AFM, 
canted, and FM states take place at n <C 1. In this case all occupied electron states lie close to the bottom of each of 
the two bands, and the momentum dependence of the electron energy can be approximated by: 



-zt-i 



2rriA 



where 



t± (the lattice constant is put to 1) and z = 6 is the number of nearest neighbors. The densities n± 



are related to the corresponding Fermi momenta pf± by 
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and the chemical potential of conducting electrons is: 
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The second part the Eq.(|7|) only holds if both electron bands are occupied. 
The energy (per site) can then be written as follows: 
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where the first term is the energy of the AFM exchange between localized spins and other terms represent the energy 
of conducting electrons (a = |(67r 2 ) 2 / 3 ). 

The large energy of the on-site FM exchange, —Jh§ti, is substracted from Eq.(||), as it only renormalizes the 
chemical potental. For n « 1, the fourth and fifth terms in Eq.(||) (the kinetic energy of electrons counted from the 
bottom of the corresponding band) are small compared to, respectively, the second and the third terms. Therefore, 
in this section the last two terms will be neglected (their role is considered in Sec. III). 

We first consider the case when the canting angle 9 is very close to tt and 7 <C 1, so that, according to Eq.(^), the 
difference between t+ and t_ is relatively small and the two electron bands are filled almost equally: 
Then, using Eqs.([I])-([7|), the energy of the QC state is approximately: 
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The optimal value of 7 can be found from the condition ^2 = an d is given by: 
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(the last expression is only valid for (n — ni)/ni <C 1). 

While in the classical treatment of de Gennes the canting of sublattices takes place at arbitrarily small doping || , 
in the quantum treatment the system stays antiferromagnetic up to the critical concentration ri\ H. This happens 



3 



because in the two-band picture electrons can propagate even in the colinear AFM state via the process described 
above, while in the classical treatment they are localized. The two bands have the same width, -^===, and are equally 
occupied, so that the energy of the AFM state is: 

Eafm = , — — JS 2 . (12) 

For n > n%, the energy of the two-band QC state becomes lower than the energy of the AFM state: 

E QC = E AFM - l J{2 Z s \ l)2 (n - n,) 2 . (13) 

As 7 grows with increasing n, so does the difference between t+ and i_. As a result, above certain concentration 
n 2 the second part of Eq.(Q) can no longer be satisfied. Thus, for n > n% only the "+" band is filled, while the "— " 
band is empty (n+ = n and n_ = 0). 

The energy of the one-band QC state is given by 

E QC = -zt+n - ~ JS 2 + zJ{2S + 1) sinh 2 7. (14) 



An optimal value of 7 for n > ri2 satisfies 
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The concentration 712, at which the second band disappears, is found using Eq.(|l5|) and 
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When 7 becomes much larger than 1, the hopping amplitde t + « t cos |, i.e., it coincides with the classical expression 
used by de Gennes. In this case we recover the CC state with the energy 

9 z 

Ecc = -ztn cos- + -JS 2 cos 9. (17) 
The optimal canting angle is given by cos | = 2 j T ^ , so that the optimal energy of the CC state (for S 1) is 

The cross-over from the one-band QC state to the CC state occurs at 

8JS 3 / 2 

n~n 3 = . (19) 



Finally, for 



MS 2 

n> n 4 = — - — , (20) 



the angle 6 — 0, i.e., the canted state transforms into the FM state with the energy: 

E FM = ~ztn+ ^JS 2 . (21) 

Summarizing the results of these section, we can say that on the class of the homogeneous trial ground states with 
two sublattices, the AFM state is the most favorable for n < n\. For rii < n < n 2 the two-band QC state has lower 
energy. For n > n 2 only one of the bands is filled. For n > n 3 the one-band QC state coincides with the CC state. 
Finally, for n > 714 the ground state is ferromagnetic. 
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III. INSTABILITY OF CANTED STATE AGAINST PHASE SEPARATION 



In this section we study the stability of the canted state against a phase separation by calculating the electronic 
compressibility, 

_ d 2 E 
dn 2 

It is easy to see that under the condition (||) and for small n the compressibility is negative. Thus, for the two-band 
QC state at n > n\ 
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[cf. Eq.(J13|)], while in the classical canted regime [see Eq.(17)] 
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Thus we see, that the homogeneous canted (both classical and quantum) state is unstable for small n. This 
instability is similar to the instability against a phase separation in the Hubbard model jl| and t-J model 

At sufficiently large n the kinetic electron energy makes the compressibility of the canted state positive. For 
arbitrary band-filling Eq. ( |l7|) for the energy of the CC state has to be substituted by 

E cc {n) = -tcos 6 -f{n)-t C os6, (24) 

where f{n) is the kinetic energy of electrons with the hopping amplitude t = 1 (which depends only on the electron 
density n), and J = zJS 2 . Then the canting angle is given by 

cos - = =- 25) 

2 2J V ; 



E CC {n)=-tf{n) 2 -i (26) 



and the energy of the CC state is 



For small band-filling f(ri) « — zn and the last equation agrees with Eq. (|18|) . While for small n the compressibility is 
negative, the increase of the electron density results in the increase of the Fermi-pressure of the electron gas, which 
stabilizes the canted state. At half-filling the electron energy, f(n), reaches minimum and, according to Eq.(p6|), so 
does the energy of the CC state. Thus, at least close to half-filling k > 0. For D = 3, the compressibility of the CC 
state is positive for n >~ 0.16. 



IV. PHASE DIAGRAM 



In this section we discuss the phase diagram of the double-exchange model. We treat spins classically (S — ► oo), 
keeping, however, J = zJS 2 finite (the numerical calculations in this section were performed for D — 3, but the 
behavior in other dimensions is qualitatively the same). In Fig. 1 we plot the energy of a homogeneous state, E(n), 
which is the CC state, for n < 77,4, and the FM state for n > n^. The energy of the CC state is given by Eq. (p6[) , 
while the energy of the FM state for arbitrary density is 

E FM {n) = tf(n) + ~. 

For small n, the CC state separates on the pure AFM phase and a second phase that contains all the electrons. 
The density of electrons in the second phase, n*, is defined by the Maxwell construction: 
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For n* > 7i4, the second phase is ferromagnetic (see Fig. la), while for < 714 it is the stable canted state (see 
Fig. lb). Which of the two situations is realized depends on the ratio of j. For the small ratio (j <^ 0.425 in 3D), 
the phase separation goes into the pure AFM and FM states, for larger j, the ground state is the mixture of the 
AFM and stable canted states. This is illustrated in Fig. 2, where we plot the phase diagram of the system. The thick 
line shows the dependence of rt* on j and separates stable states from unstable states. In the ferromagnetic region 
of the diagram grows with j, while in the canted region w 0.291 independent of j (both sides of Eq.(|^) with 

E = Ecc can be divided by y , after which the equation becomes independent of t and J) . In the unstable FM region 
and a part of the unstable canted region the compressibility is positive, so that the corresponding states are locally 
stable. They are, however, globally unstable against the phase separation, as follows form the Maxwell construction. 

In the J}j t and S > 1 limit, we consider in this section, the phase diagram of the double-exchange model is 
symmetric around n = -|. This can be seen by applying to the Hamiltonian (Q) the particle- hole transformation 




+h\ for i £ A 
-h\ for i £ B 



where ct annihilates electron on site i with spin parallel to the local spin on the same site, while hi is a corresponding 
hole operator. This transformation leaves the kinetic energy term (first term in (Q)) unchanged, but it reverses the sign 
of the Hund's coupling (the second term in (Q)). The fact that holes are coupled to local spins antiferromagnetically, 
rather than ferromagnetically, does not, however, make any difference: the important thing is that the orientation of 
the spin of the hole is uniquely determined by the orientation of the local spin. Thus, the double-exchange system 
with n holes per site (e.g., with electron density 1 — n) has the same properties as the system with n electrons per 
site. 

Experimentally, however, there is a marked difference in the behavior of the hole-doped manganties (Lai_ x M x Mn03, 
M = Sr, Ca, x < 0.5) and the electron-doped (overdoped) manganites (x > 0.5, the electron density related to doping 
by n = 1 — x). While the properties of the hope-doped systems qualitatively resemble those described above, the 
electron-doped manganites (with predominantly Mn 4+ ) behave quite differently: they are mostly insulating, often 
with charge ordering in the form of stripes fill and sometimes with anisotropic magnetic structures |l6|] . The reason 
for the assymetry of the phase diagram is not completely clear. One of the factors, which is not taken into account 
in the double-exchange model studied above, but which may play an important role, is the orbital degeneracy of the 
electrons bands p7| . Note, however, that the tendency to a phase separation discussed in this paper is also present 
for the degenerate bands |l7j . 



V. POLARON STATE 

There exist many different realizations of an inhomogeneous phase-separated state. For the model (|l|) the most 
favorable one is the complete phase separation (considered in the previous section), in which case one part of the 
sample contains all the charge carriers and is ferromagnetic (or strongly canted), while the other part is an undoped 
antiferromagnet. Such a phase separation would, however, be counteracted by the Coulomb interaction (not included 
in the Hamiltonian (|l|)). The Coulomb interaction may suppress the phase separation completely, or it could favor 
the formation of the FM droplets (spin-bags) immersed into the AFM matrix. More complicated structures, e.g., 
stripes |l8| are also possible. In general, the phase separation in the double-exchange model in the presence of 
Coulomb interactions seems to be rather similar to that in the t-J model considered in the context of the high-T c 
superconductors |^JT^|. In the later case, the instability towards the phase separation remains even with the Coulomb 
interaction taken into account. 

In this section we consider a simple type of the phase separated state: the AFM matrix with the FM "droplets" , 
each containing one electron. In other words, we assume that electrons propagating in the AFM background form 
ferromagnetic polarons (HJ7|- Such a state, with only one electron per droplet and with the droplets well separated 
from each other will help to minimize the Coulomb energy contribution mentioned above. We will demonstrate that 
the energy of this state is, indeed, lower than that of the homogeneous canted state. 

We calculate now the energy of the system for small doping, assuming that each electron is localized inside 

a FM sphere of radius R. The energy of such a state can is given by 

Ep = - - tn (z - + n~nR 3 - - f 1 - u-ttrA . (28) 
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The first term in ( p§| ) is the energy of the lowest state of an electron in the sphere, while the second and the third 
terms give, respectively, the exchange energy of localized spins in the ferromagnetic regions and in the remaining 
antiferromagnetic matrix. The ferromagnetic polarons are assumed to have sharp boundaries and the surface energy 
is neglected. Furthermore, we assume that the polarons do not overlap, in which case the change of the energy of the 
system due to doping is simply proportional to n. 

The optimal polaron radius, found from the condition ^ = 0, is given by 

-or 

and the energy of the polaron state is: 

E P = -l- ztn + ^n(7rf) 3/5 (2J) 2 / 5 . (30) 

It is easy to see that the energy of such a polaron state is lower than the energy of both the classical and the 
quantum homogeneous canted states. Moreover, at n — the polaron state is also more favorable than the FM 
state. Only at higher concentration, 
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the transition to the saturated FM state occurs. At this concentration polarons begin to overlap and for small j this 
concentration up to a numerical factor coincides with discussed in the previous section. 

Ferromagnetic polarons can be also formed in materials with layered electronic and magnetic structure, such as 
manganites. In LaMnC>3 the layers are formed by planes of ferromagnetically coupled localized spins, while the 
exchange between the spins from two neighboring layers is antiferromagnetic. Also, the interlayer and intralyer 
hopping amplitudes, tu and tj_, arc different. Such a layered magnetic and electronic structure is a result of the 
orbital ordering that takes place in this material JlSj] . The polaron is elongated ("cucumber") for t\\ > t±_ and 
compressed ("pancake") for tn < t±. A detailed discussion of the FM polaron in the layered system will be given 
elsewhere. 



VI. CONCLUSIONS 



In conclusion, we have shown that the tendency to a phase separation and a formation of a spatially inhomogeneous 
state, found for several models of strongly correlated electrons, is also an inherent property of the double-exchange 
model. For large Hund's coupling, Jh > ( > JS 2 , and small doping the homogeneous canted state has negative 
compressibility. This is related to the instability of the canted state against the phase separation into pure antiferro- 
magnetic and ferromagnetic phases. The latter phase contains all charge carriers. For zJS 2 ~ t the compressibility 
of the canted state is negative only at small concentrations of charge carriers. In that case, it can separate into the 
pure antiferromagnetic phase and the canted state with a large concentration of charge carriers and a large canting 
angle. 

The spatially inhomogeneous state may occur in a variety of different forms, e.g., a random mixture of AFM and 
FM phases, magnetic polarons ("fog" of small FM droplets), or some regular structure, such as stripes. Which of this 
situations is realized would depend on specific conditions and may be different in different parts of the phase diagram. 
However, the tendency to form an inhomogeneous state seems to be a generic property of the double-exchange model. 
The specific features of manganaties (layered magnetic structure, orbital degeneracy), though definitely important for 
determining a particular form of the ground state, nevertheless seem to preserve this tendency. 

Several theoretical problems still remain to be clarified. The first one, already mentioned above, is the role of the 
Coulomb interaction in determining a detailed nature of the inhomogeneous state (the average size of polarons and 
their spatial ordering (cf. experimental results p5[)). Another question is the detailed spin structure both inside 
and outside polarons. In the treatment above we made a simplifying assumption that the inner part is the saturated 
ferromagnet, while the outer part is the pure antiferromagnet without canting. However, in general, this need not to 
be the case: both these phases may have certain degree of canting. In particular, since the spectrum of spin excitations 
in AF matrix is gapless, the distortion of the perfect antiferromagnetic order should decay slowly with the distance 
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from a polaron || . All these questions definitely require further study, but at this point we can, in any case, say that 
the homogeneous canted state seems to be nearly always unstable, and inhomogeneous structures of some kind have 
to be realized in low-doped double exchange systems, including CMR manganites. 

Turning to an experimental situation in manganites, we can state that there exist already many indications that 
the phase separation is indeed present there. The formation of stripe phases in overdoped Lai-^Ca^MnOa (x > 0.5) 
p5| can be viewd as a manifestation of this tendency. However, much more indications of this are now observed 
in the most interesting region of x < 0.5 and, in particular, in the CMR regime. Already the old data |2l]] on a 
coexistence of both the FM and AFM phases in the magnetic neutron scattering, which are sometimes treated in a 
picture of homogeneous canted state, are actually much better interpreted as a mixture of FM and AFM microregions 
P|. Recent neutron scattering results p0| confirm this picture and, moreover, show the presence of spatial liquid-like 
correlations of the FM droplets (in the system studied in p(| these were not purely FM droplets, but rather the 
microregions with the large canting angle). Direct confirmation of the existence of the inhomogeneous state comes 
from the NMR measurements |]l2| , which show the presence of two different hyperfine fields corresponding to FM and 
AFM regions. And finally, quite recently the presence of an inhomogeneous state with a random mixture of phases 
with different electronic properties, was directly visualized by the STM study of Lao.7Cao.3MnC>3 [^2|: the metallic 
and insulating regions with an average scale of 10 — 50nm were directly seen in this investigation. 

The natural consequence of the phase separation is that the main properties of CMR manganites, including transport 
properties, should be treated in the percolation picture |ic|| . This concept explains quite naturally the critical doping 
x c ~ 0.16, above which the ferromagnetic metallic state appears in Lai-^M^+MnOs |jjo[ . An experimental support for 
the percolation picture comes from the recent study which shows that I— V characteristics of Lai-^Pr^Cao.aMnOa 
are strongly nonlinear even at small voltages: such behavior is expected for the percolated state, for which the electric 
field would modify the distribution and shape of metallic clusters and enhance conductivity. 

Thus, both from the theoretical treatment and from many experimental results it follows that, most probably, the 
electronic and magnetic state of doped manganites is intrinsically inhomogeneous, leading to the percolated state. 
This should be definitely the case for small doping (for which the random potential of impurities should be also taken 
into account); but it is also quite probable that such is the nature of even optimally doped manganites (x ~ 0.3) with 
CMR. This gradually emerging picture, promisses to give a new explanation of the properties of these fascinating 
materials. 
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Figure Captions 



Fig 1. The energy of a homogeneous state £ as a function of electron density n (thick line) calculated for t = 1 
and two different values of J: ~ 0.344, corresponding to n 4 = 0.2, (a) and <~ .484, corresponding to n 4 = 0.4. The 
Maxwell construction (thin line) gives the density of the phase containing all the electrons in the phase-separated 
state. For n 4 < n*, the phase-separated state is the mixture of the AFM and FM states (a), while for ra 4 > it is 
the mixture of the AFM and CC states (b). 

Fig 2. The phase diagram of the double-exchange model. The states to the left of the thick line are unstable 
against the phase separation. The thick line shows the dependence of the density of the electron-reach phase on 
J. The thin line separates the canted state from the ferromagnetic state. The regions of the CC state with positive 
and negative compressibility are separated by the dotted line. 
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